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^ Abstract 

O, We use AdS/CFT to construct the gravitational dual of a 5d CFT in the back- 

<^ ground of a non-extremal rotating black hole. Our boundary conditions are such 

that the vacuum state of the dual CFT corresponds to the Unruh state. We ex- 
tract the expectation value of the stress tensor of the dual CFT using holographic 
' _(' renormalisation and show that it is stationary and regular on both the future and 

■*— ' the past event horizons. The energy density of the CFT is found to be negative 

Q^ everywhere in our domain and we argue that this can be understood as a vacuum 

Zh polarisation effect. We construct the solutions by numerically solving the elliptic 

I— I Einstein-DeTurck equation for stationary Lorentzian spacetimes with Killing hori- 

^^ zons. 
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1 Introduction 

Arguably, the formulation of the AdS/CFT correspondence [H [21 13] has been the most 
important advance in theoretical physics in recent years since it provides a definition of 
quantum gravity in asymptotically anti-de Sitter (AdS) spacetimes in terms of the dual 
conformal field theory (CFT). Therefore, this correspondence provides a mean to address 
fundamental questions in quantum gravity in a well-defined setting. One of the most 
remarkable effects in quantum gravity is Hawking's realisation that black holes evaporate 
|1]. Hawking's result motivated the study of quantum fields in black hole backgrounds 
(see e.g., [3] for a detailed discussion) but in most of these calculations the fields are 
assumed to be non-interacting. Adding interactions is very difficult in this approach. On 
the other hand, in AdS/CFT we can study strongly interacting quantum fields in curved 
black hole backgrounds by constructing classical solutions to the bulk Einstein equations 
with a negative cosmological constant. 

The authors of El E] have pioneered the study of Hawking radiation in large N^ 
strongly-coupled field theories which admit gravity duals. As these references pointed out, 
the vacuum state in which the CFT is defined is specified by the boundary conditions that 
the geometry satisfies deep in the IR. Then, in principle, the AdS/CFT dictionary allows 
one to extract all physical quantities of the CFT in that particular black hole background 
in the chosen vacuum state. 

Using the numerical methods developed in [9], Ref. ^U\ constructed an Einstein 
metric in 5c? with a negative cosmological constant such that the boundary geometry 
is the 4d asymptotically fiat Schwarzschild solution. In the IR the geometry reduces 
to the Poincare horizon of AdSs and therefore, by the AdS/CFT correspondence, this 
gravitational solution should describe the large N^, strong coupling limit of A/" = 4 super 
Yang-Mills (SYM) on the background of the Schwarzschild black hole in the Unruh state 
[TO]F] This solution allows one to extract the leading 0{N^) piece of the expectation value 
of the quantum stress tensor of the dual field theory in this particular vacuum state. The 
results of [lO] showed that, remarkably, (Tj^) is static, regular on both the future and 
the past event horizons, and there is no energy fiux at infinity. These results contradict 
the free field theory expectations (see e.g., [12]), showing that strong interactions can 
significantly alter the behaviour of quantum fields, especially in curved spaces. Using 
the construction of [lO], Ref. [13j subsequently addressed a very closely related problem, 
namely the construction of a static black hole localised on the brane in the single infinite 
Randall-Sundrum II braneworld model O US]. In this paper, static braneworld black 
holes of both small and large radii relative to the radius of the parent AdS space were 
constructed. These results contradicted the non-existence conjecture of [161 113 tlH] , which 
is based on free field theory intuitionjj 

The description of CFT's in black hole backgrounds in the Hartle-Hawking state in 
terms of their gravity duals requires different boundary conditions in the IR. In the Hartle- 
Hawking state, the stress tensor of the plasma at infinity should approach that of a thermal 
fiuid in equilibrium with the boundary black hole. From the bulk point of view, deep in 
IR the geometry should asymptote to the planar black hole. Ref. [6j further predicted 
the existence of two families of solutions depending on whether the bulk horizon has 
a single connected component (black funnels) or two disconnected components (black 



^This reference did not distinguish between the Unruh and Boulware vacua but [llj has pointed out 
that the Boulware state requires a minimum energy configuration and hence an extremal horizon. 
^Refs. [ini [20] had remarked that free field theory intuition might not hold. 



droplets). The black funnels describe a situation in which the plasma couples strongly to 
the boundary black hole, resulting in an efficient exchange of heat between them. On the 
other hand, the black droplets describe the weakly coupled situation. These two phases 
should be connected through a Gregory-Laflamme type transition. The black funnels have 
been subsequently constructed in [TT] (and their properties conffimed) whilst the black 
droplets remain elusive so far (see [6l [71 El IH] for more details). 

Black funnels have motivated further studies of heat flow in AdS/CFT. Ref. [21] 
has shown that rotating BTZ black holes in AdSa can be reinterpreted, by means of a 
change of conformal frame, as describing heat flow in the CFT. However, the conformal 
symmetry in two dimensions implies that the left and right temperatures of the CFT 
must be constants, and the bulk horizons are Killing. On the other hand, [22] more 
recently considered a stationary funnel solution in global AdS4 consisting of a single bulk 
horizon connecting two boundary black holes. The latter have different temperatures 
which implies that there is a net transport of heat along the bulk horizon, which cannot 
be Killing]^ 

In this paper we generalise the solutions of [TU] by adding rotation to the boundary 
black hole. Whilst in principle the Hartle-Hawking state cannot exist in rotating black 
hole backgrounds with Killing horizons due to issues with superradiance [12] , we are not 
aware of any obstructions for the Unruh state. There is an extensive literature on studying 
free quantum fields in static black hole backgrounds (see [S] and references therein) but 
we are not aware of any results in rotating black hole backgrounds in four or higher 
dimensions. In this paper we remedy this and by using AdS/CFT we compute, for the 
ffist time, the expectation value of the stress tensor of a strongly interacting 5d quantum 
field theory (and in particular a CFT) in a rotating black hole background. Our boundary 
conditions are such that the vacuum state of the CFT is the Unruh vacuum. 

For technical reasons that will be explained subsequently, we construct the gravita- 
tional dual of a 5d CFT in the background of the 5d Myers-Perry black hole with equal 
angular momenta [2^ . However, our construction should capture most of the qualitative 
physics of the Unruh state in rotating black hole backgrounds for strongly coupled CFT's 
that admit gravity duals. Another motivation for our work comes form the recent inter- 
est in studying 5d super conformal field theories (SCFT's). In general, five dimensional 
gauge theories are non-renormalisable and therefore do not exist as microscopic theories. 
However, there are a few notable exceptions. Indeed, for certain choices of gauge group 
and matter content there exist certain supersymmetric gauge theories in 5d that arise 
as superconformal fixed points of the renormalisation group [251 ISHl |27]. One example 
of these gauge theories that is relevant for this paper is 5d super Yang-Mills with gauge 
group USp(2A^c) and matter consisting of a hypermultiplet in the antisymmetric repre- 
sentation of the gauge group and Nj hypermultiplets in the fundamental representation. 
This theory arises as the low energy worldvolume theory on Nc D4 branes probing an 08 
plane with A^^ D8 branes on top of it [25] . The theory admits a supergravity dual [271 l2H] 
in lOd type IIA massive supergravity. Ref. [22] computed the free energy of this field 
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theory on a S^ using localisation techniques and showed that it scales like Nc . Moreover, 
the field theory calculation agrees with the supergravity result, which uses the calculation 
of entanglement entropy [30] , generalised to gravity duals with a non-trivial dilaton [31] . 
More recently, Refs. [321 E3] generalised these theories to quiver-type gauge theories and 
studied their supergravity duals. 



^Black holes with non-Kilhng horizons describing stationary plasma quenches were constructed in |23) . 



The rest of this paper is organised as follows. In ^ we explain our setup and provide 
the details of the numerical construction of our solutions. In ^ we present our results. 



In ^3.1 we consider some geometrical aspects of our solutions and the holographic stress 



tensor is discussed in §3.2[ Section ^ is devoted to the conclusions and directions for 
future research. Further details of our solutions as well as convergence tests are relegated 
to Appendices |A] and |B] respectively. 

Note added: While our work was nearing completion, we became aware of [3l], which 
has appeared simultaneously with our paper on the arXiv. |3l] also addresses the con- 
struction of the gravitational dual of bd CFT's in the background of the bd Myers-Perry 
black hole with equal angular momenta. 

2 Rotating boundary black holes in Poincare AdS 

In this section we describe our construction of rotating black holes in Poincare AdS. More 
precisely, we will consider Einstein metrics in Qd with a negative cosmological constant 
such that the boundary metric is conformal to the bd Myers-Perry black hole ^24^ with 
equal angular momenta. In the region of the spacetime which is far from the horizon, 
our geometry approaches the Poincare horizon of AdSe. Therefore, by the gauge gravity 
duality, our solution is dual to a CFT on a rotating black hole background. As discussed 
in [Tot [TT] , our boundary conditions in the IR should imply that the dual CFT is in the 
Unruh vacuum state. 

Since the boundary horizon is rotating (and Killing), by continuity the bulk horizon 
will have the same angular velocity as the boundary black hole. Therefore, our solutions 
are the first examples of rotating horizons in the Poincare patch of AdS q Moreover, we 
shall restrict ourselves to the case of non-extremal black holes, and the temperatures of 
the boundary and bulk horizons are the same. 

The reason why we consider the equal angular momenta Myers-Perry black hole as our 
boundary geometry is a technical one: such a spacetime is of co- homogeneity one, and 
therefore the full spacetime metric will depend only on two coordinates (see ^2.2 below). 



The problem then reduces to that of solving partial differential equations (PDE's) in two 
variables. We believe that our solution should capture the essential physics of quantum 
field theories on rotating black hole backgrounds. In some respects, the properties of 
the bd Myers-Perry solution with equal angular momenta resemble those of the Ad Kerr 
solution {e.g., existence of a smooth extremal limit, non-existence of an ultraspinning 
regime [35]), and therefore we expect that the physics of TV = 4 SYM on the Kerr 
black hole background should be qualitatively similar to what is described in this paper. 
Nevertheless, it would still be interesting to study the latter. In this other case, the 
problem is of co-homogeneity 3 (but still tractable with our techniques) and we will leave 
it for future work. 

We want to construct a solution to the Einstein equations in Qd with a negative 
cosmological constant obeying the boundary conditions outlined above. To do so, we 
will use the method of |i9i i36j (see ^37j| for a review) and instead of solving the Einstein 
equations, we will solve the so called Einstein-DeTurck equations (also referred as the 
'Harmonic' Einstein equations), 
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-Rab = -Rab + ^ fi'ab + V(a6) = , C = 9 " i^" be - '^°' be) ^ (2-1) 



"^The planar limit of a rotating black holes in global AdS yields a boosted black brane. 



where I denotes the radius of AdSp 



r"^^ is the Levi-Civita connection associated to the 



spacetime metric g that we want to determine. For simphcity, we take F"^^ to be the Levi- 
Civita connection associated to some chosen reference metric g on the same spacetime 
manifold M.. Solutions to (2.1) are Einstein iff ^" = everywhere on M.. To achieve this, 
we shall impose boundary conditions (see ^2.2.1) which are compatible with ^"'\dM = 
[SI Uni [36] and a posteriori we find that indeed ^" = everywhere in our domain (within 
our numerical accuracy). 



2.1 Equal angular momenta bd Myers-Perry black hole 

In this subsection, we begin with a review of the geometry of the bd Myers-Perry black 
hole with equal angular momenta. We then introduce new coordinates which will be 
useful for our subsequent numerical construction of the gravitational duals of CFTs in 
rotating black hole backgrounds. 

The bd Myers-Perry black hole with equal angular momenta is described by the fol- 
lowing line element 






dt2 + 



{R^ + a^Y - liE? 



di?2 + 



/x 



B? + a^ 



dt + -a^ 



{R' + a') dfi?3) , 
(2.2) 



where fi and a are the mass and angular momentum parameters respectively, and 

cr^ = — sin -ipdO + sin 9 cos ^d0 



0"^ = cos ipd9 + sin 9 sin t/jdcj) 
a^ = cos 9d(f) + dip 



(2.3) 



are the left-invariant one-forms on SU(2). Here d^2?3^ = | ((a^)^ + (cr^)^ + (cr^)^) is the 
round metric on the unit 3-sphere. The event horizon is located at i? = Rh where Rh is 
the largest real root of the equation {Rj^ + a^)^ — ^Rj^ = 0- ^^ the following it will be 
useful to express the mass parameter /i in terms of Rh and a: 



/^ 



{Rji + a 
Rjj 



2\2 



(2.4) 



It is then straightforward to find a convenient compact radial coordinate by defining 

2 Rh + 0. 



R' + a' 



'1 — r 



2^2 ' 



(2.5) 



so that r = corresponds to the event horizon (or, more precisely, the bifurcation surface) 
and r — )■ 1 corresponds to spacelike infinity. In terms of the new radial coordinate r, we 



can rewrite (2.2) as 



dsl.. = - r'Toir) dt' + -^^^^, dr' + — ^^^ {a' - fio(r)dt)' + 



'MPs 



4 1 



On 



4(1 



:^ ^^(2) 



(2.6) 



where 



Ro{r 
So{r 
To(r 

Boir 



2)(i?|^-a2(l-r2)2) 



Rl + a^ 






(2.7) 



.2\4\ 



/?2 



2afl-r 



2a4 



R]j + a^ 



n4' 



and dr2?2) 



cr 



1n2 



a ) is the metric on the round unit 2-sphere. Note that when 



written in this form, the metric exhibits the structure of a U(l) fibration over an 5* base. 



As we shall see in ^2.2, we can consistently construct an extension of this metric into the 



bulk of AdSe whilst preserving this structure. In these coordinates, the null generator of 
the future event horizon is given by. 



k 



d_ 
dt 



Q 



H 



d_ 

dip' 



(2i 



and the angular velocity and surface gravity of the horizon are 

2a o To(0) {R 



n 



H 



fio(O) 



Rl + a^ 



K 



4i?o(0) Rl{Rl + a^ 



(2.9) 



The form of the metric (2.6) is particularly useful because the boundary conditions at 



the horizon that the various functions satisfy become transparent (see ^2.2). In fact, any 
rotating black hole solution with a Killing horizon and with an orthogonally transitive 
isometry group can be cast in an analogous form, so the boundary conditions that we 
discuss below are general. Alternatively, we could have chosen to work in co-rotating 
coordinates as in [3S] at the expense of losing the manifest asymptotic flatness of the 
metric. This could potentially be problematic when extending the metric into the AdS 
bulk if we want to require that the spacetime asymptotes the Poincare horizon of AdSg 
in the region far from the horizon. 



2.2 Metric ansatz and boundary conditions 

Let us consider first the metric on the Poincare patch of AdSg: 



•-l^AdSe 



^(dz^ - dt^ + dR^ + R^dny^) 



(2.10) 



For our construction it is convenient to introduce new coordinates which are adapted to 
the Poincare horizon of AdS. To do so, we generalise the coordinate change introduced in 
ITOl and define new coordinates (r, x) as 



{l-x^)^R^ 



l-r2 



R 



v/(2 - x^) {Rl + a^) 



l-r2 



(2.11) 



in terms of which (2.10) becomes 



ds 






+ -. i-T + — ^ — TT^ dnj 



(2-x2)(l-r2)2 Rl + a^ {I - r^)^ {I - r^f ^^^ 

(2.12) 
In these coordinates the conformal boundary of AdS is in the asymptotic region x — ?• 1, 
and X = is the fixed point set of the SO (4) symmetry. These coordinates are adapted 
to the Poincare horizon of AdS in the sense that the latter hes at r = 1 and x < 1. The 
point X — !■ 1 and r — )■ 1 corresponds to spacehke infinity at the boundary. 

We want to construct a 6-dimensional asymptotically locally AdS spacetime such that 
its boundary metric is in the same conformal class as the 5d Myers-Perry black hole, 
and that approaches the Poincare horizon of AdS far from the conformal boundary. The 
isometry group of the boundary metric is MiXSU(2)xU(l). It has been shown in 4c? 
that regular and static asymptotically locally AdS metrics inherit the isometry group 
of the boundary metric [M]- We shall make the reasonable assumption that in higher 
dimensions and in the presence of rotation the same result holds true. Therefore, our full 
spacetime metric should have the full ]RixSU(2)xU(l) isometry group of the boundary 
metric. Then, the most general metric {i.e., closed under diffeomorphisms that preserve 
the symmetries) that satisfies our symmetry assumptions is given by 



ds 



2 



'1 — r 



2^2 



4:X(x,r) , , 2xr(R% + a'^)Y(x,r) , , 



{1-x^Y L(2-x2)(l-r2)" (l-r2 

-r To{r)T{x,r}dt H — —^ [a^ - ll{x,r)dt) 



x'^{2-x^)Bo{r)B{x,r) 
4(l-r2)2 
x^ (2 - x^) So{r) S{x, r) . 4 Ro{r) R{x, r) 






(2.13) 



where {X{x,r),Y{x,r),T{x,r),B{x,r),Q{x,r),R{x,r), S{x,r)} are the unknown func- 
tions to be solved for, and {To(r), i?o(r), i?o(?"), 5'o(r)} are the Myers-Perry functions 



given in (2.7). As above, the coordinate ranges are < x < 1 and < r < 1 but 
we will effectively treat the asymptotic regions x = 1 (conformal boundary of AdS) and 
r = 1 (Poincare horizon) as boundaries and work on the square domain < x < 1 and 
< r < 1. On this domain we shall assume that all our unknown functions are smooth 
and bounded. Note that because we have chosen coordinates which make the isometries 
of the metric manifest (and hence reducing the computational cost), r = (event hori- 
zon) and X = (fixed point set of the SO (3) symmetry) appear as boundaries of our 
domain. In fact, we can identify the square {0 < r < 1, < x < 1} with the orbit space 
O = A^/[]R(XSU(2)xU(l)], which has the natural structure of a manifold with bound- 
aries and corners. However, as discussed in [TOl |36], some of these boundaries are in fact 
only fictitious boundaries after suitable smoothness conditions on the metric functions 
are imposed there. 

2.2.1 Boundary conditions 

Refs. [ini EH] discuss in generality the behaviour of the metric on different kinds of 
boundaries that are relevant for our problem. In the following we shall specialise these 
general results to our case. 



Comparing (2.13) with (2.6) we see that at x = 1 we have to impose the following 



Dirichlet boundary conditions: 

X{l,r)=f. 



F(l,r) = 0, fi(l,r)=fio(r) 



T(l,r) = i?(l,r) = 5(1, r) = 5(1, r) 



This choice ensures that as x — > 1 the metric becomes 



(2.14) 



R]j + a? 



ds' 



[l-r 



2n2 



'l-x^ 



(2-x2 



dx' 



{Rjj + a 



ds2 



(2.15) 



which shows that the boundary metric lies in the same conformal class as the 5d Myers- 
Perry black hole. It is worth mentioning that we allow for the constant factor of (i?^ + a^) 
multiplying ds^p^ so that the boundary conditions at x = 1 and r = 1 (see below) are 
compatible for any choice of Rh and a. This is necessary to ensure smoothness of the 
unknown functions. 

Following the general discussion on boundary conditions for extremal horizons in [10] , 
since at r — )■ 1 we want our metric to smoothly approach the Poincare horizon of AdSg, 
we must impose 



X 

R 
Q- 



+ (l-r)Xi 

e2 



Rl + a^ 



, Y = {l-r)Y,, T 
r)Ri, B 



;i-r)fii, Ti\r=i-Ri 



r=l 



Rl + a^ 
= constant 



Rl + a^ 
+ (l-r)Bi 



+ {l-r)T^ 



S 



R2 



+ (l-r) Si, 
(2.16) 



where all the functions with the subscript i are smooth at r 
that i?o(l) = Bo{l) = 5*0(1) = Rjj + a^ and also that dQ^^-^ = 



1 and X < 1. Note 



dQf^^ + (a 



3^2 



As 



shown in [TU], these boundary conditions are necessary and sufficient for our spacetime 
to asymptote the Poincare horizon of AdSg. In fact, we shall not impose the condition 
Ti|r=i — -Ri|r=i = constant but we check that it is satisfied on our solutions as consequence 
of the equations of motion. Hence, we can use this condition as an estimate of the 
numerical error of our solutions. 

At X = both the S"^ of the base space and the U(l) fiber degenerate. To ensure that 
they do so smoothly, the functions X, Y, T, R, B, Q, S must be smooth functions of x^ near 
X = and hence they should obey a Neumann boundary condition there. In addition, to 
avoid the presence of conical singularities at x = we must require 



So(r)5(0,r) So{r)S{0,r) 



X(0,r) 



X(0,r) 



1 at X = 



(2.17) 



Smoothness of the manifold at the horizon requires all metric functions to be smooth in 
r^ near r = 0, and hence they obey a Neumann boundary condition there. In addition, 
as discussed in |36], we must ensure that our metric has the same surface gravity and 
angular velocity as the reference metric (see below for our specific choice of reference 
metric), which can be achieved imposing 



T(x, 0) = i?(x, 0) , fi(x,0)=fio(0) 



(2.18) 



Then the surface gravity and angular velocity of our spacetime will be given by (2.9). 



In order to solve the Einstein-DeTurck equations (2.1) we have to specify a reference 



metric g on the same manifold M. as the physical spacetime metric g that we are seeking. 
This is equivalent to choosing a particular gauge. In our case, we have chosen g to be 



given by (2.13) with 



X 



r = 0, 



R = S 



Rl + a? ' 



^H 



\L = 1 in , B 



Rl + a^ 



x^ + 



RUl 



X 



Bo{Rjj + a^ 



(2.19) 



Finally we note that the boundary conditions spelled out above are compatible with 
^""IdM = 0. In the static case and for boundary conditions similar to ours, Ref. [TU] 



proved (under suitable smoothness assumptions) that all solutions to (2.1) obeying these 
boundary conditions should satisfy ^" = everywhere on Ai and therefore are solutions 
to the Einstein equations as well. An analogous result that applies to stationary metrics 
is not available yet so, given our boundary conditions, we cannot rule out a priori the 
existence of 'Ricci solitons' {i. 
the problem is elliptic, the solutions should be locally unique 



solutions to (2.1) with ^" ^ 0). However, because 

Therefore an Einstein 



metric cannot be arbitrarily close (in the space of solutions) to a Ricci soliton P, [101 [36] . 
Therefore, we can always check a posteriori for a given solution whether ^" = 0. Needless 
to say, for the solutions presented in ^we find ^" — )■ in the continuum limit. 



2.3 Details of the numerics 



In this subsection we present the details of our numerical solution of (2.1) with the 
boundary conditions specified in the previous subsection. 



To solve (2.1) we considered the associated parabolic equation, namely the Ricci- 



DeTurck flow equation, and evolved the metric according to 



d_ 



9ab 



"2 Rab 5 



(2.20) 



until we reached a fixed point and hence also a solution of (2.1). Here A is the Ricci flow 
time. In order to carry out the flows we have to provide an initial metric gab\x=o and for 
the results presented in ^we chose the initial metric to be equal to the reference metric. 
The results should not be qualitatively different for other choices of initial data, but we 
have not investigated this. 

As discussed in |39l |9] , the convergence of the flow towards the fixed point depends on 
the stability of the fixed point. The latter is determined by the presence of negative modes 
in the spectrum of the Lichnerowicz operator computed around the fixed point jj In our 
simulations we found that our initial data always converged to the desired Einstein metric 
with no fine-tuningjj indicating that our solutions are stable fixed points of the Ricci- 
DeTurck flow and hence do not posses negative modes. Even though the 5d Myers-Perry 
black hole has a negative mode |35j (completely analogous to the celebrated negative mode 



The linearisation of (2.1 1 around an Einstein metric coincides with the Lichnerowicz operator so in 



this discussion we are assuming that the fixed points are Einstein metrics. This is indeed the case for the 
solutions presented in this article. 

^See [2] for a detailed explanation of how to use Ricci flow as a method to construct Einstein metrics 
possessing one negative mode. 



of Ad Schwarzschild |10] or Kerr [H]), it seems that pinning the metric to the boundary 
of AdS projects out this negative mode. Ref. [10] observed the same phenomenon for 
the 5d bulk case in which the boundary metric is simply 4d Schwarzschild. It would be 
interesting to have a more detailed understanding of why and how this happens. 

In AdS/CFT only the conformal class of boundary metrics matters and therefore we 
can choose any representative within the class without loss of generality. For the problem 
at hand, choosing a representative of the conformal class amounts to fixing the overall 
length scale of the boundary 5d Myers-Perry black hole. In the results presented below, we 
fixed the mass parameter n to some convenient value and varied the angular momentum 
parameter a to move along the family of solutions. Note that this implies that both the 
surface gravity and angular velocity vary along this family. 

As in [9], there are two main sources of numerical error. The first one is present in 



any numerical method and it is simply the fact that we discretise the equations (2.20) 
according to some scheme and seek a solution to the discrete problem. The latter ap- 
proximates the continuum solution and, for a consistent discretisation, in the continuum 
limit the error should decrease at some rate which depends on the discretisation scheme 
and the degree of differentiability of the solution. The second source of error is related 



to the fact that it takes infinite Ricci flow time for (2.20) to reach its fixed point, but 



in practice we can only evolve (2.20) for a finite time. We then have to make sure that 
we are sufficiently close to the fixed point so that the dominant error is the discretisation 
error. We do so by monitoring |£^i?/30 + l|max (which vanishes for an Einstein metric in 
the continuum limit) along the fiow and we only stop the simulation when this quantity 
is constant within some tolerance. The value of this constant gives an estimate of the 
numerical error of our solution and it decreases as we increase the spatial resolution; in 
Appendix [B] we present more details about the convergence tests that we have performed. 
Of course, the amount of Ricci flow time that it takes to reach the flxed point depends on 
how close to extremality the solution is and the resolution of the spatial discretisation. 

We simulated the Ricci-DeTurck flow using the method of lines. The two spatial 
dimensions were discretised onto a square grid of Chebyshev points, and spatial derivatives 
were approximated by spectral differentiation. The flow was then time-evolved using a 
third-order Runge-Kutta integrator with a flxed step size. We note also that spatial 
differentiation is equivalent to simple pointwise multiplications in the spectral space, and 
therefore we can improve the computational speed by using the discrete cosine transform 
rather than repeated matrix multiplications. 

In our calculation of the curvature tensor, we hard-coded the expression for each of 
the non-zero components in terms of the metric components and their derivatives. Due to 
the complexity of this problem, a typical component involves a large number of additions 
and subtractions of floating point quantities. As the metric becomes singular near the 
horizon, the relative sizes of these summands can become vastly different. At double- 
precision, we found that this becomes problematic when we get sufficiently close to the 
fixed point: numerical error starts to dominate and the fiow oscillates wildly and does 
not converge any further. The problem worsens as we increase the number of grid points 
since the metric components need to be computed closer to the horizon. 

The problem is exacerbated by the fact that, in order to extract the boundary stress 
tensor, we need to read off the fifth derivative from our numerical data, and for each extra 
derivative that we require, the numerical results needs to be evaluated more precisely. 
In some cases, we found that in order to extract the fifth derivative to an acceptable 
quality, we had to invoke quadruple-precision fioating point arithmetic in our numerical 



computation. Unfortunately, as no current equipment implements this at hardware level, 
we had to rely on software libraries to emulate higher-precision floating point arithmetic, 
which can be over 20 times slower than native computation. 

To optimise the time that it takes to obtain solutions with acceptable numerical accu- 
racy, we divided the typical run into three phases. We first evolve the Ricci flow equation 
numerically using the standard 64-bit double-precision floating point arithmetic. Once 
our error indicators (Ricci scalar and norm of DeTurck vector) start to oscillate, we took 
the numerical solution at the final Ricci flow time and use it as initial data for a second 
run, which is carried out using 80-bit "long double" precision. Finally, we take our solu- 
tion from the second run as the initial data for a final run at 128-bit quadruple precision. 
Typically it takes between one and two days to complete the whole process. In all cases, 
the data presented in the next section is the result of this final run at quadruple precision. 

For the majority of cases, we found that a 36 x 36 square Chebyshev grid provided 
a good balance between accuracy and computational requirements. However, for the 
slowly-spinning cases we also found that finer grids are required in order to maintain an 
acceptable accuracy (we chose a 45 x 45 grid for the three most slowly- spinning cases). 

3 Results 

In this section we present our results: §3.1 discusses some aspects of the geometry of the 



solutions and in §3.2| we present the holographic stress tensor together with its interpreta- 
tion. To construct our solutions we have fixed the mass parameter of the boundary black 
hole, fi, and then varied the angular momentum parameter a such that the dimensionless 
surface gravity is K/ia = ^ with n = 1, . . . , 16. Therefore, for our solutions the horizon 
of the black hole is always non-extremal. 

3.1 Geometry of the horizon 



To visualise the horizon geometry of our solutions {i.e., fixed points of (2.20)) as we vary 
the angular momentum of the boundary black hole, we consider the following embedding. 
We take the geometry of the spatial cross-sections of the horizon transverse to the U(l) 
fiber, 



ds' 



H n _ rr2\2 



(1-X2 



4X(M) 2 , x^(2-x^)go(0)g(2;,0) ^^, 
2_a;2 ^^ + ^ ^^^(2) 



(3.1) 



and embed them into Ad hyperbolic space with the same radius £, ds^^ = ^{dz^ + di?^ + 
E^dVt^,^))^ ^s a curve R = R{z). Then we colour-code this curve according to the relative 
size of the U(l) fiber with respect to the base S"^, which is measured by the function 
a{z) = Bo{0)B{x,0)/{So{0)S{x,0)) where the coordinate x should be regarded as being 
a function of the coordinate z of the ambient hyperbolic space. The freedom in the 
embedding is fixed by requiring that for the static case the radius of the boundary black 
hole is -R(O) = 1, which fixes the ADM mass of the boundary black hole for any a > 0. 
The results are depicted in Fig. [l] As this figure shows, the radius of the S'^ becomes 
smaller as angular momentum increases; this is simply the reflection of the fact that the 
horizon area at fixed mass decreases as the angular momentum increases. Also, note that 
the extent of the horizon into the bulk decreases as the angular momentum increases. 
Further, for a given fixed value of the angular momentum parameter a > 0, the squashing 
of the horizon three-sphere is more severe at the boundary rather than at the interior of 
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Figure 1: Embeddings of the radius R{z) of the horizon 5*^ into four- dimensional hyper- 
bohc space for fixed mass parameter fi of the boundary black hole. At the boundary of 
hyperbolic space, z = 0, the value of this radius function coincides with the corresponding 
one in the 5d Myers-Perry black hole geometry. The colour code depicts the squashing of 
the horizon S^ according to the relative size of the fibre with the respect to the base S"^, 
which is measured by the function a{z) = Bo{0)B{z,0)/{So{0)S{z,0)). The blue curve 
with R{0) = 1 corresponds to the zero rotation case and a increases moving inwards in 
this plot. The squashing becomes more severe near the boundary and with increasing 
angular momentum of the boundary black hole. Absence of conical singularities at the 
axis of symmetry R{z^s.x) = implies that there is no squashing there. 



the spacetime. In fact, the regularity condition (2.17) at the axis of symmetry, x = 



implies that the the horizon S^ must be locally round there, and therefore there cannot 
be any squashing. 

In Fig. [2] we depict the evolution of the curvature invariant CabcdC"'''^'^ i'^ along the flow 
for a typical solution. This is a useful geometric quantity because it gives information of 
both the curvature of the spacetime and the correctness of our boundary conditions. As 
this figure shows, Cabcd.C'^^'^'^ l^ never blows up along the flow (and in particular at the 
fixed point), which suggests that there are no curvature singularities in our domain. In 
addition, we see that for any A this curvature invariant vanishes at both the Poincare 
horizon of AdSe and at the conformal boundary. This confirms that our spacetime indeed 
has the correct asymptotics near these two boundaries. 



3.2 Stress tensor 

We extract the vev of the stress tensor of the dual CFT from the bulk gravity solution 
using the standard holographic renormalisation prescription [12]. For a Qd bulk spacetime, 
the expectation value of the stress of tensor the dual CFT is given by 



{T.,) 



levrGf 



(5) 



(3.2) 
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Figure 2: Cabcd.C"'^'^''' l'^ along the Ricci flow in the whole domain for the K/ia = 0.15 
solution. The snapshots correspond to the Ricci flow time A = 0, 0.5, 1 and the fixed 
point (from top to bottom in this plot). After Ricci flow time of A = 1 the surface 
becomes almost indistinguishable from that of the actual fixed point. As this plot shows, 
Cabcd.C"''"^'^ l^ does not blow up anywhere in the domain along the flow. Also, it vanishes at 
both the Poincare horizon and at the boundary of AdS, which indicates that our spacetime 
is asymptotically AdS near these boundaries. 

(5) 

where gij denotes the 5th order term in the near boundary expansion of the metric in 
Feffer man- Graham coordinates and Gq is the Qd bulk Newton's constant. 

Our bulk Qd solutions can be trivially uplifted on half an 5"^ (see [13]) to IIA massive 
supergravity in lOrf, which arises as the low energy limit of intersections of D4 and D8 
branes in Type I' string theory. The CFT that arises as the fixed point of this system 
under RG flow is USp(2A^c) Sd SYM with matter consisting of a hypermultiplet in the 
antisymmetric representation of the gauge group and Nf hypermultiplets in the funda- 
mental [28]. This is the CFT dual to our supergravity solutions. As we have reviewed in 

5 

ql , the number of degrees of freedom of this theory scales like Nc at large N^. and hence 



(3.2) should scale with the same power of Nc. One way to see that this is indeed the 



case is the following. Consider the string frame metric of our solution in lOd massive IIA 
supergravity [43j|^ 



^■5(10) 



sm a 3 



- [eds%) + ^(da2 + cos2 a ds% 



(3.3) 



where ds?gj is our Qd asymptotically AdS solution and a G (0, vr/2], so that only the upper 



^The lOd solution has a curvature singularity that is not hidden behind the horizon. However, for 



-3/10 



large N^ there is large a region of the spacetime with sin a^ Nc 
dilaton are small and thus one can trust the supergravity solutions 



where both the curvature and the 
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hemisphere of a full S^ is covered. Of the other supergravity fields we will only need the 
dilaton|3 

= ^^^ (sin «) ' ■ (3-4) 

Ref. [22] showed that the relation between the radius of AdSe, i, and the string length 
for the class of theories that we are considering is 

n- 8-Nf- ^'-'^ 

and the lOd Newton's constant is taken to be 2kIq = {2TTy g'^i^^ In the decoupling limit 
the low energy excitations are sensitive to the string frame metric and therefore it seems 
natural to relate the lOd Newton's constant to the 6d one by dimensional reduction on 
the S"^ in the string frame. By doing so we take into account the fact that the dilaton is 
non-trivial. This gives the following relation between the lOd Newton constant and the 
6d one: 



1 1 3(8 -Nf)l^/N: 



2nl 2 Kfo 272' 



TT 



/ da A vol^s (sinQ;)3 cos^ a . (3.6) 



Combining the result of this integral with (3.5) gives 



£^ 27 



leyrGe A<dy/2T^^ ^8 - N^ 



(3.7) 



5 
vc • 



which shows that indeed (3.2) scales like Nc 

Now we proceed to extract the expectation value of the stress tensor of the dual field 
theory using holographic renormalisation. Since we are interested in spacetimes with 
]RtxSU(2)xU(l) isometry group, we can consistently truncate the general form of the 
line element in Fefferman-Graham (FG) coordinates to this particular class of metrics. 
Solving the Einstein equations in the usual near boundary expansion we find that the 
stress tensor is given by 

(T,,) dxMx^' = -^^ [T^{R) de + R^{R) dR^ + B,{R){a' - n,{R) dt)^ + S,{R) dfiJ^J , 

(3.8) 
where T5, R^, B^, Q^, S5 can be extracted from our numerical solutions. These functions 
are not all independent: the Einstein equations at 5th order in the near boundary expan- 
sion (in FG coordinates) impose an algebraic constraint between them which is nothing 
but the tracelessness condition of the dual stress tensorFj At 6th order one finds a dif- 



ferential constraint which implies that (3.8) is covariantly conserved. Therefore, we find 



that the stress tensor is fully specified in terms of three independent functions. However, 
in our construction we shall independently extract the five functions from our numerical 
data, then use the trace and divergence of (Tij) to estimate the associated numerical error. 
We find that the trace is largest near the horizon, with magnitude of around 10~^ to 10^^ 



*In field theory, a necessary condition for a fixed point to exist is Nf < 7 |55]. In supergravity, tfie 
solution describing the D-brane system only makes sense if this condition is satisfied [28] . 

^Recall that 2k^ = leirGd- 

^"Recall that in odd number of boundary dimensions, the stress tensor is always traceless regardless of 
the boundary metric being Ricci flat or not because there is no gravitational conformal anomaly. 
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(the fast rotating cases giving smaller values). This quickly decays to well below 10~^ in 
all cases as we move away from the horizon. 

In order to obtain the components of the stress tensor in terms of our numerical 



solutions we have to change from the working coordinates (x, r) used in our ansatz (2.13) 



into the Fefferman-Graham ones, {z,R) that we used to get (3.8). We proceed as in 
and determine the change of coordinates z = z{x,r), R = R{x,r) in a near boundary 
expansion, requiring that ^" = at each order in (1 — x^). Note that in our setup ^" has 
only two non-vanishing components so imposing .^'* = order by order determines the 
coordinate change completely. 

The dual stress tensor (T^) is symmetric in its two indices but because the boundary 
metric (2.2) is not positive definite, the linear map (T*) from vectors to vectors need 



not be diagonalisable. However, in our case we find that it is, and for all values of the 
angular momentum parameter a of the boundary 5d Myers-Perry black hole that we have 
considered there is one timelike eigenvector and four spacelike eignvectors. We can identify 
the corresponding eigenvalues with the energy density and the pressures in the following 

way. We write the stress tensor as 

{T)) = p{R) fmj +Pi{R) [{s'y ® {s')j + {s^y ® (s^),] +P2{R) {s^ ® {s^)j +p3{R) R'®Rj 

(3.9) 
where f is the unique timelike eigenvector normalised so that t*tj = —1, and (s'')* and 
i?* are the orthonormal spacelike vectors in the hd Myers-Perry backgrouncp] written 



in Boyer-Lindquist-type coordinates, (2.2). In this way we identify p{R) as the energy 



density of the plasma seen by a local observer with velocity t* and the pj(i?)'s are the 
corresponding pressures,^ 

In Fig. |3] we plot the energy density for different values of the rotation parameter of 
the background. First note that the energy density is negative everywhere in our domain 
for all a > 0. Moreover, the energy density is a monotonically increasing function of the 
radial coordinate {i.e., it becomes less negative for large R) and it is largest (in absolute 
value) at the horizon. Further, as the angular momentum increases (moving up in this 
figure) the energy density becomes more negative at the horizon, but it is always finite. 
By continuity the energy density should also be finite at extremality. Next we consider 
the pressures (see Fig. [6] in Appendix |A] for the plots). As in |23] we find that the 
radial pressure is positive whilst the pressures along the angular directions are negative, 
and therefore there is an anisotropy in the system. Note that this anisotropy is already 
present in the static limit. Adding rotation to the boundary black hole breaks the full 
S0(4) spherical symmetry of the static case down to SU(2)xU(l), which results in some 
further anisotropy between the fibre and the sphere directions. This simply reflects the 
squashing of the S^ of the boundary black hole. As angular momentum of the boundary 
black hole increases, the horizon S'^ becomes more squashed which gives rise to a more 
negative pressure along the fibre direction than along the base 5*^ directions. 

In Fig. |4] we study the behaviour of the energy density as a function of the radial 
coordinate R in the i? — )■ cxd region. As one can see from this figure, the energy density 
decays like l/i?^ near infinity and the other components of the stress tensor exhibit the 
same asymptotic behaviour (see Fig. [T] in Appendix lAl). Not surprisingly, this can be 
understood as in the Ad case [H]. The energy density stored in the field around a point 



^^We choose (s^)^ ^ {o'^)i with the obvious proportionahty factors and Ri = \/ Trt+o^W^^JIW i'^^)i- 
^^Recall that if fc* is a unit timehke and future directed vector field, then — T*- fc-' represents the energy 
density seen by an observer with velocity given by A:°. So —p{R) is the eigenvalue of the matrix T* . 
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Figure 3: Plot of (minus) the energy density for different values of the rotation parameter 
versus the Boyer-Lindquist radial coordinate R. We have used a logarithmic scale on the 
y axis in order to enhance the near horizon region. As this plot shows, the energy density 
is negative everywhere and it becomes more negative as the rotation of the black hole 
increases. 

5 

particle in 5-spacetime dimensions behaves like oc Nc Rg/R'^, where Rg is the gravitational 
radius of the point particle. The pressures are comparable in magnitude to the energy 
density and they also decay hke 1/i?^. It should be possible to reproduce this result 
by adapting the linearised calculation of [45J to the present situation, where the bulk 
spacetime is 6-dimensional. 

Finally we consider the rotation of the dual plasma. To do so, we write the unique 
timelike eigenvector of the stress energy tensor of the dual CFT as 



at 



n{R) 



d_ 

dip 



(3.10) 



This allows us to identify VL{R) as the angular velocity of the plasma with respect to a 
static observer at infinity in the boundary directions. In Fig. [5] we plot this quantity 
for different values of the angular momentum parameter a. As one might have expected, 
VL{R) — )■ VIh as i? — )■ Rh- That is, at the horizon the plasma is co-rotating with the 
black hole. This is an example of frame dragging even though at the boundary gravity is 
non-dynamical. This had to be the case since the boundary black hole is rotating rigidly 
with respect to infinity; as there is no flux of CFT plasma through the horizon then it 
must be that the plasma is co-rotating with the black hole. In the asymptotic region, 
R — )■ oo, we find that the angular velocity of the plasma decays like l/R"^, so that at 
infinity the plasma is at rest and the timelike eigenvector coincides with the asymptotic 
timelike Killing vector. 
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Figure 4: Plot of the rescaled energy density for K/ia = This plot shows that the 

energy density decays like 1/W in the asymptotic -R — )■ oo region. The other components 
of the stress tensor have a similar behaviour near infinity. The colour code is the same as 
in Fig|3] 

We now discuss the physical interpretation of our results. From Fig. |3] we see that 
the energy density of the plasma is negative. Because the dual stress tensor is traceless, 



the weak and the strong energy conditions are equivalent in this case and therefore (3.9) 



violates the usual energy conditions considered in classical general relativity. The same is 
true for the stress energy tensor of A/" = 4 SYM in the background of Ad Schwarzschild in 
the Unruh vacuum computed in [10] . However, the energy density in the Hartle-Hawking 
state computed in [TT] turns out to be positive. One may be concerned that the fact that 
we are seeing a negative energy density may be a sign of a pathology and/or an instability 
of these solutions. However, this need not be the case since our stress tensor corresponds 
to that of a quantum field theory and as such it need not satisfy any of the classical energy 
conditions. 

Ref. [TU] (see also [H]) provided a heuristic physical interpretation of the stress tensor 
(3.9) as a 'halo' of plasma in equilibrium with the boundary black hole. According to 



this picture, the strong attractive self-interactions of the CFT are balanced by thermal 
radiation pressure from the black hole. There is another instance in which putting a 
QFT in a non-trivial background results in a stress tensor with some properties that 
are qualitatively similar to the ones we find. This is the case of the celebrated Casimir 
effect [12] , in which the vacuum energy of the electromagnetic field between two perfectly 
conducting plates is negative, giving rise to an attractive force between the two plates. 



Therefore, the fact that the energy density in (3.9) is negative may simply reflect the 



attractive nature of the self-interactions of the CFT and the non-vanishing vev of the 
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Figure 5: Rescaled angular velocity of the plasma as function of the radial coordinate R 

for Kfi2 = Near infinity the angular velocity of the plasma decays as ~ 1/R'^ and at 

the horizon it reduces to the angular velocity of the black hole. 



stress tensor may be due to the Casimir effect 
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We may now try to better understand (3.9) in terms of the Casimir effect or equiva 



lently vacuum polarisation. To do so we may concentrate on the t-r part of the geometry 
since this gives rise to the non-trivial causal structure of the spacetime: the sphere direc- 
tions are mere spectators in both this paper and in [10] because of symmetry assumptions. 
Therefore, we may try to understand the behaviour of the t-t and r-r components of the 
stress tensor by considering a simple 2d model. For concreteness we will consider (T*) 
for a generic quantum field theory in the 2d space 



ds' 



1 - 



2M 



dt^ + 



dr^ 



2M ' 



(3.11) 



where we have taken the t-r part of the 4d Schwarzschild solution but our considerations 
should apply without change to the t-r part of our geometry. The advantage of working 
in 2d is that the stress tensor is completely determined by symmetries and conservation 
in terms of purely geometrical data |17l HE] . Indeed, assuming that the stress tensor is 
time-symmetric, one finds that the conservation equation reduces to 



driK) 



M 



r^ 1 — 



2M 



{{T\) - (n,)) 



dr 



1-^1 (n> 



M 

^{T\)- (3.12) 



^^Ref. [T7] noted a large contribution of the Casimir effect in the stress tensor of a conformally coupled 
weakly interacting scalar field in the background of the BTZ black hole. 
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This equation can be readily integrated to get 



(T.>=diag|-i% + (r,) ^^('' 



1 _ 2_M ' \ «/' 1 _ 2M 
r r 



(3.13) 



where the first entry denotes the t-t component of the stress tensor and 

j{T\){r') 



Hor) 



M 



dr 



2M 



„/2 



(3.14) 



Therefore, we see that the full stress tensor is completely determined by its trace. The 



latter is fixed by the conformal anomaly, which in 2d and for (3.11) gives 

AaM 



(3.15) 



where a is a constant that depends on the spin of the field being considered. For a free 
massless scalar field one has a = ^^ JITJIIS]- Using this in (3.14) allows us to completely 
determine the expectation value of the full stress tensor in this geometry: 

(n>="d-8(^-I^(l + ¥)(l + *J?).I^(l + ^)(l + ^))- (3-16) 



In the near horizon region (3.16) has the same qualitative features as our 5d stress 



tensor or that of [lOj. In particular, it is static and regular on both the future and the past 
event horizons, and the energy density near the horizon is negative. Clearly, the origin 



of (3.16) is due to the spacetime curvature near the horizon and hence one may conclude 



that it is a vacuum polarisation effect. Of course, for large r (3.16) reduces to the stress 



tensor of radiation at the Hawking temperature set by the horizon in (3.11), and clearly 



this is not what we see. Therefore, this 2d model captures only the physics in the near 
horizon region but not in the asymptotic region near infinity. 

Now consider adding the extra dimensions to take into account the full Ad Schwarzrschild 
geometry. In this case, we shall assume that the stress tensor is static, spherically sym- 
metric and traceless because we are interested in considering CFTs in this background. 
With these assumptions, the conservation of the stress tensor fully determines it up to 
one function [IHl ID] : 



(T)) = diag 



Hjr) 



2M ^ 
r ' 



2e(r), 



BaW) 



HAr) 



r2(l 
dr'(r'-3M)0(r'), 



r2(l 



2M \ ' 



e(r), 0(r) 



(3.17) 



2A/ 



where the first two entries in (3.17) correspond to the t-t and r-r components and 
(T^g) = (T'^) = 0(r) by spherical symmetry. Adding the extra dimensions introduces 
new pressures along the angular directions and the extra powers of r in the denominators 



in (3.17) can alter the behaviour of the components of the stress tensor near infinity. 



The function 0(r) in (3.17) determines the full stress tensor and its actual form de- 



pends on details of the CFT under consideration. However, one can put some constraints 
on it on general grounds, without assuming strong coupling. First, note that in order for 
(3.17) be finite at the horizon (r = 2 M), we only need to require that 0|r=2M be finite. In 
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addition, we have seen above that the 2d model captures the vacuum polarisation effect 
near the horizon and therefore we need to require that (T*^) > for r k.2 M. This implies 
that < in this region and hence we are led to conclude that the polarisation of the 
vacuum near the horizon will inevitably make the pressures along the angular directions 



negative. Also note that from (3.17) one concludes that at the horizon r = 2 M the follow- 
ing relations must hold: {T\)\r=2M = {T^r)\r=2M = — 0|t-=2A/- A close examination of the 
data of [in] shows that this is true for A/" = 4 SYM in the background of Schwarzschild, 
the deviations being smaller than 1%. At infinity we want (T*) — )■ and for this we need 



B ~ 1/r or faster. Then (3.17) implies that all the components of the stress tensor have 
the same behaviour near infinity. Note that B has to decay sufficiently fast if we want 
that {T\) and {T\) do not change sign at a sufficiently large r. One can infer the fall 
off of B from the linearised calculation of [35] and one finds that the pressures along the 
sphere directions are negative and the fall off is ~ 1/r^. 

Finally we note that even though we find negative energy densities, quantum energy 
inequalities (see [32] for a review with references) can be used to bound them. Quantum 
energy inequalities are the remnants of the classical energy conditions satisfied by quantum 
field theories and they quantify the extent to which a quantum field can violate these 
classical energy conditions. Even though most of the results available in the literature 
apply to Minkowski spaces, on physical grounds one expects that similar results hold on 
curved spacetimes over regions which are small compared the curvature length scales. In 
fact, recent progress [SO] shows that this intuition is indeed correct. For instance, consider 
a globally hyperbolic spacetime Ai with a timelike Killing field t"' admitting a foliation 
into constant time surfaces such that A4 ~ M( x S. For any region Eq C S such that 
the spacetime metric takes the Minkowski form, the energy density of a free scalar field 
satisfies [5'OJ^ 

(T.,n'^nX>-(^^, (3.18) 

for any Hadamard state u, and where n" is the unit vector along t". Here r{x) is the 
radius of the largest Euclidean 3-ball which can be isometrically embedded in Sq and 
C4 = 3.169858..., is a constant. For Dirac or vector fields one can obtain similar type of 
inequalities. Our results, obtained using AdS/CFT and valid at strong coupling, show 
that in this regime the energy density also remains bounded from below. 

4 Conclusions 

In this paper we have constructed the gravitational duals of 5d CFTs in rotating black 
hole backgrounds. For simplicity we have considered the case in which the boundary black 
hole is the 5d Myers-Perry solution with equal angular momenta because the metric is of 
co-homogeneity one, and hence the full spacetime is of co-homogeneity two. Moreover, in 
the region far from the horizon the spacetime reduces to the Poincare horizon of AdS and 
hence the CFT should be in the Unruh state. Whilst our solutions should capture the 
general qualitative physics of the Unruh vacuum in rotating black hole backgrounds (at 
large Nc and strong coupling) it would still be interesting to consider the gravitational 
dual of A/" = 4 SYM on the background of Kerr. From a computational viewpoint the 
problem reduces to solving elliptic PDEs in 3 variables and we hope to address it in 



^■^ Interactions can also be included and they do not change in essence the quantum energy inequalities, 
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the future. The construction of the gravitational duals of CFTs on rotating black hole 
backgrounds becomes even more interesting if one insists on considering an IR horizon 
with finite temperature. As p!Tj pointed out, this requires the horizon of the black hole 
be non-Killing and recently methods for constructing this kind of solutions have been 
developed [2SII22]. 

We have also noted that a generic feature of the Unruh vacuum for strongly coupled 
CFT's in black hole backgrounds is that the energy density is negative everywhere in the 
domain of outer communications. We have argued that this can be understood in terms 
of the polarisation of the vacuum due to the curvature of the background. Further, our 



stress tensor (3.9) is bounded from below, in accordance generic local quantum energy 
inequalities j;49j. The fact that the energy density is negative need not signal that our 
solutions are unstable but we have not performed a detailed study of their stability under 
gravitational perturbations. It would be very interesting to carry out such a study. 

It is interesting to compare the results of this paper (and ^lOj) with the black funnels 
[TT] . The latter represent the gravitational dual of the Hartle-Hawking state for A/" = 4 
SYM on the background of Schwarzschild. In the black funnels case, the energy density 
is positive everywhere and not only at infinity (see Fig. 4 in [H]), where the stress tensor 
reduces to that of pure radiation at the Hawking temperature set by the IR horizon. 
Therefore, their stress tensor does not seem to have a component corresponding to the 
vacuum polarisation. This should be contrasted with the free field theory results fSiT] , 
which indicate that near the horizon it is the vacuum polarisation that dominates the 
stress tensor, as in this paper. However, [6] conjectured that there should exist another 
family of solutions which are dual to the Hartle-Hawking state, namely the black droplets. 
The latter should arise as a continuous deformation of the solutions of [10], where the 
continuous parameter is the temperature of the IR horizon. Because in [TU] (and of 
course in our present solutions as well) the non-vanishing stress tensor is due to vacuum 
polarisation effects, the stress tensor of the black droplets should be dominated by vacuum 
polarisation effects near the horizon, at least when the temperature of the IR black hole 
is sufficiently small compared to the size of the boundary black hole. In the asymptotic 
region, the stress tensor should reduce to pure radiation, as in the black funnels case. 
Therefore, for the black droplets the {T\) component of the stress tensor should change 
sign, being positive near that horizon and negative in the near infinity region. It would 
be very interesting to construct these solutions. 

Acknowledgements 

We are grateful to Sungjay Lee, Nick Dorey, Diego Rodriguez- Gomez and especially Har- 
vey Reall and Toby Wiseman for discussions. PF supported by an EPSRC postdoctoral 
fellowship [EP/H027106/1]. ST is supported by the European Research Council grant no. 
ERC-2011-StG 279363-HiDGR. The computations presented in this paper were done on 
COSMOS at the UK National Cosmology Supercomputing Centre in Cambridge. 

A Further details of the stress tensor 

In this appendix we collect more results about the two remaining eigenvalues, namely 
Pi{R) and P2{R), which together with the energy density p{R) discussed in the main text, 
fully determine the stress tensor. 
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In Fig. [6] we plot the absolute value of the two pressures, although both of them are 
negative everywhere in our domain (see Fig. [?]). In these plots we have used a log scale 
in the y axis so as to zoom into the region near the horizon. As shown in this figure, both 
pressures are finite at the horizon. It is interesting to note that the pressure along fibre 
direction is always larger (in absolute value) than the pressure along the S"^. In Fig. u^we 
plot the rescaled pressures. As this plot shows, both pressures decay as ~ 1/ W at large 
distances, as advertised in the main text. 

B Convergence tests 

When a numerical scheme is employed to solve a continuous problem, it is essential to 
verify that that the solution so obtained converges appropriately as the continuum limit 
is approached. For spectral methods, the error should decay exponentially as the number 
of grid points is increased, provided that the function being approximated is smooth. In 
this appendix, we present evidence to show that our numerical construction does indeed 
exhibit this behaviour. 

We seek a metric which satisfies the condition Rab = — 5 gab-, where we have set the AdS 
radius to one, C. = 1. By taking the trace of this equation, we see that the scalar curvature 
should approach R = —30 in the continuum limit. Therefore, we use the maximum of the 
quantity |i| + l| as a measure of error in our solution. In Fig. 8 we plot this quantity for 
the metric which is the fixed point of Ricci flow over a range of grid sizes. We performed 
this convergence test for two parameter values K/iS = 0.0625, 0.75. In both cases, we can 
see that the error follows a general exponentially decaying trend, even though the decay 
is not monotonic in the grid size. 

Furthermore, our a posteriori gauge-fixing scheme allows for the possibility that the 
flow converges towards a Ricci soliton. To check that our solution is indeed Einstein, we 
also track the decay of the squared-norm of the DeTurck vector: = C,"'C,a- Again, we 
expect that this quantity should become zero in the continuum limit. Fig. [8] shows that 
|0| also decays exponentially as we increase the number of grid points. 

We should also comment on the behaviour of these quantities as we evolve the Ricci 
flow. In Fig. M we plot the maximum of the scalar I ^ + 1 1 for a typical evolution as a 
function of Ricci flow time. Here we can clearly see the flow stabilising as we approach the 
fixed point. The figure also shows how insufficient precision causes the flow to oscillate 
before the fixed point is reached. In all cases, continuing the flow using higher-precision 
floating point arithmetic allows the flow to progress further. The corresponding plot for 
|<^|max exhibits the same behaviour. 
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0.05 




Figure 6: Pressures along the sphere directions (top) and along the fibre {bottom) in the 
near horizon region. Notice the log scale of the y axis. The pressures are finite and 
negative everywhere in the domain. 
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Figure 7: Rescaled pressures along the sphere directions (top) and along the fibre (bottom). 
Both pressures decay as 1/R'^ for _R — )• oo. 
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Fierure 8: LoErarithmic plot of I ^ + ll and Icbl as a function of the grid size. The red 

^ ^ -^ I oU I max ' ^ ' max ^ 

and purple curves correspond to /t/iS = 0.0625 and 0.75 respectively. In all cases the error 
decays exponentially, but the more highly-spinning case has a larger decay coefficient. 
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Figure 9: Logarithmic plot of U^ + 1 as a function of Ricci flow time for a typical 
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run of the simulation. The simulation starts in the blue portion at double precision, then 
continued using long double precision in the red portion, and finished off using quadruple 
precision in the green portion. 
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